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STRAIN-RATE EQUATIONS FOR CALCULATION OF SECONDARY CREEP 
DEFORMATION OF THICK-WALLED TUBES WITH 
INTERNAL OR EXTERNAL PRESSURE 
by Richard E. Morris 
Lewis Research Center 

SUMMARY 

The design of tubular components of nuclear propulsion systems for long service 
life at high temperature and pressure requires that the secondary creep deformation be 
controlled or limited. Examples of such structures include heat exchangers and 
nuclear -reactor fuel-pin cladding. 

Heat exchangers must operate with internal pressure at high temperature for long 
lifetimes. An important requirement is the limitation of the creep expansion of the 
thick-walled heat-exchanger tubing. 

In the case of nuclear-reactor fuel-pin cladding, both external coolant pressure 
and the development of internal fission gas pressure cause creep deformations that 
must be limited. 

Strain-rate equations are presented for use in these design problems. The stress 
and strain-rate relations apply to thick-walled tubes for both internal and external 
pressure. Both temperature and pressure were assumed to be constant. Stress and 
strain rate were assumed to be related by a power law which includes a function of 
temperature as a parameter. Thus, the equations are applicable to the stress and 
strain-rate behavior of tubes over a range of temperatures. 


INTRODUCTION 

Tubular components of nuclear propulsion systems may be designed to operate 
under extreme conditions of temperature and pressure such that creep deformation 
occurs continuously during operation. 



An important limitation in designing heat exchangers for high-temperature opera- 
tion is the creep strength of heat-exchanger materials. Similarly, the creep strength 
of nuclear -reactor fuel-pin cladding materials poses a limitation on the operating tem- 
perature of the fuel pins in a reactor . 

Helium -to-air heat exchangers may operate with an internal helium pressure of 
7 to 14 meganewtons per square meter at temperatures from 1000 to 1150 K for life- 
times of 50 000 hours. Under these conditions the heat -exchanger tubes undergo con- 
tinuous creep deformation. Rupture of heat- exchanger tubes must be avoided. The 
amount of deformation accumulated during the operating lifetime of the tubing must be 
limited by design. 

Tubular nuclear -reactor fuel- pin cladding contains the nuclear fuel. Fission gases 
are generated as the fuel is consumed. Gases collect in the high-temperature core 
void of a fuel pin. The gas is constrained by the volume of void provided in the fuel- 
pin design so that pressure develops as a function of fuel burnup. Pressures as high as 
25 meganewtons per square meter may be generated during the life of a pin. 

For high-performance reactors the peak temperature of fuel-pin cladding may be 
from 1100 to 1500 K. In this temperature range, the regularly increasing fission gas 
pressure will cause the cladding to creep during the 10 000-hour design life that is 
typical of these fuel pins. Radioactive fission gases must be contained within the fuel 
pin. Rupture must be avoided. Total creep deformation during the lifetime of the fuel- 
pin cladding must be limited by design. 

The coolant in the gas-cooled reactor must operate under pressure sufficient to 
obtain the heat-transfer characteristics needed to cool the fuel pins. For helium, a 
coolant pressure of about 10 meganewtons per square meter is needed. This pressure 
is sufficient to cause continuous compressive creep of the fuel-pin cladding during the 
first part of the life of the reactor. Compressive creep must be limited to avoid col- 
lapse of fuel-pin cladding. The fission gas pressure gradually increases to oppose the 
coolant pressure. Fission gas pressure is greater than the coolant pressure at the end 
of the life of the reactor. 

Each of these design problems involves the relation between stress and strain rate 
in thick- walled tubes undergoing creep from either external or internal pressure. 

This report presents stress equations and strain-rate equations for applications to 
these problems. The analysis follows the method presented by Bailey (ref. 1) which is 
a three-dimensional stress analysis of internally pressurized tubes undergoing second- 
ary creep. The analysis presented in this report changes one of his fundamental 
assumptions so that the strain-rate equations apply to a range of operating tempera- 
tures. In addition, by using the boundary conditions for external pressure the analysis 
yields stress equations and strain-rate equations for tubes under external pressure. 

All of the equations apply only to the case of constant pressure and temperature. 
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SYMBOLS 


A material constant, hr - ^(N/m^) -11 

a inside radius, cm 

1 O n 

B material constant, hr (N/m ) 

b outside radius, cm 

C constant of integration 

D constant of integration 

AH apparent activation energy, J/(K)(mole) 
n stress exponent 

9 

p pressure, N/m 

R gas constant, 8. 3143 J/(K)(mole) 

r radius, cm 

T absolute temperature, K 

u radial displacement, cm 

e creep rate, hr 

€ equivalent creep rate, hr~* 

P b/a 

2 

ct stress, N/m 

2 

ct equivalent stress, N/m 

X b/r 

Subscripts: 

i internal 

o external 

r radial 

z axial 

0 circumferential 
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ANALYSIS 


This analysis follows the work of Bailey for the three-dimensional stress analysis 
of internally pressurized tubes (ref. 1) except for the power law assumption that re- 
lates stress and strain rate. Bailey assumed a power law e = Ao* 1 . In this form both 
A and n are temperature dependent. Consequently, strain-rate relations for creep 
in thick-walled tubes under internal pressure based on this power law are applicable 
only for a relatively narrow range of temperature for which A and n are constant. 

The power law assumed in this analysis contains a temperature correction term 
-AH/RT 

e ' Use of this term reduces the temperature dependence of the constants in 
the stress, creep-strain-rate relation. With the temperature correction term included, 
the creep-strain-rate relation used in this analysis is 

i = BcrV iH / RT 


This power law has been used by several investigators, for example, Maag and Mattson 
(ref. 2). Garofalo (ref. 3) indicates that this strain-rate relation applies to secondary 
creep at constant temperature for many metals and alloys at low stresses. Long life- 
times with low creep rates at high temperature limit the use of many metals to low 
stresses. 

The use of the power law with the temperature correction term in this analysis and 
adding the boundary conditions for tubes with external pressure provide strain-rate 
relations for thick-walled tubes loaded with static internal or external pressure under- 
going isothermal secondary creep as a function of temperatures. 

The assumptions used in the analysis are as follows: 

(1) The tube material is isotropic. 

(2) Creep associated with the deformation of the tubes is secondary creep. 

(3) Secondary creep strain follows the von Mises flow rule. 

(4) Stress and strain rate are assumed to be related by the empirical power law 

<f = ±Ba n e -AH//RT . (The sign is selected for tensile (+) or compressive (-) strain rate. ) 

(5) The axial creep rate is assumed to be zero. Planes perpendicular to the axis 
of the tube remain plane. The creep deformation is one of plane strain. 

(6) The. principal axes of stress and creep strain rate are coincident and remain 
so during Strain up to moderate amounts. 

(7) The principal shear strain rates are proportional to the principal shear stresses. 

(8) Poisson’s ratio is assumed to be 1/2. This assumption requires that the sum 

of principal orthogonal creep rates equal zero, that is e r + + e z - 0. 
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INTERNAL PRESSURE 


The use of the power law assumption with the temperature correction factor results 
in strain-rate equations for and e r . Correlation of the equations with experi- 
mental data will provide values for the material constants. The resulting equations 
are more general than Bailey’s strain-rate equations, in that the equations presented 
here include temperature as a parameter. 

The stress equations derived by Bailey (ref. 1) are not affected by this change in 
the stress-strain rate power law assumption. The equations for the stresses in a tube 
undergoing secondary creep at constant temperature are given here in the terminology 
of this report for completeness (the derivation is included in the appendix): 


a -^ p . 

r p 2 /” - 1 1 


1 - 


M x 2/n 


a e=' 


n 


pVn _ ! 1 


1 - SjJ\ vVn 


a_ = 


n 


p 2 / n - 1 


- l/3 x 2/n 

a . ^ Pi 

n p 2/n - 1 


(1) 

(2) 

(3) 

(4) 


Equation (4) is obtained by substituting the stress equations into the relation for the 
maximum distortion energy stress (equivalent stress) 


a = -± 




K - a e )2+ (<T e - a z) 2 + K ~ a r )2 


1/2 


Use of the modified power law assumption yields the following strain-rate equations: 

L = - (5) 

r 2 
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(6) 


e „ = l^BaV AH / RT 

6 2 

e = 0 (assumption (5)) (7) 

e = Ba n e~^^' // ^^' (assumption (4), tensile creep rate) (8) 

These eight equations define the stresses and strain rates in an internally pressur- 
ized tube at constant temperature for secondary creep conditions. 


External Pressure 

The derivation for stresses in tubes from internal pressure uses the boundary 
conditions 


cr = -p at r = a 

n r = 0 at r = b 

Changing the boundary conditions to 

CT r = 0 at r = a 

a r = -p at r = b 

we have the condition for tubes loaded by external pressure. Then if the boundary con- 
ditions are used to evaluate the constants of integration in the derivation, the equations 
for the stresses become 


2/n _ 2/n 


cr = - ^ ^ p 

r n /_ F o 


p 2/n - 1 


a 9 ~ 


p 2/n . x 2/n 


p 2/n - 1 


(9) 


( 10 ) 
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( 11 ) 


p2/n _ (2^1U 2/n 


a z = - 


n 


p 2 / n - 1 


--V5- x 2/n p 

n p 2 / n - 1 


( 12 ) 


Assumption (4) is modified for negative strain rates assuming that compressive 


behavior is similar to tensile behavior. The resulting strain rates are 

e = Be”'^ H//RT o^ 1 (13) 

r 2 

£ -_^3 B e- AH / RT a n ( 14 ) 

9 2 

e = 0 (assumption (5)) (15) 

z 

7 = -Ba n e” AR / RT (assumption (4), compressive creep rate) (16) 


Equations (13) to (15) relate the stresses and strain rates to functions of pressure, 
geometry, temperature, and materials properties for a given radius in the wall of a 
tube. The sign of the empirical exponential stress function a n is assumed to be posi- 
tive. 


DISCUSSION 

Equations are presented for the stresses and strain rates in isothermal thick- 
walled tubes with static internal or external pressure that are applicable over a range 
of temperatures. Tube structures that are designed for long lifetimes at high tempera- 
tures may undergo creep deformation continuously. Design stresses must be selected 
so that creep rates are limited to very low secondary creep rates. If this is done, ter- 
tiary creep deformation will be avoided. The equations presented apply only for tubes 
undergoing secondary creep deformation. 

The power law between stress and strain rate used in developing the equations con- 
tains a temperature correction term. Thus, the equations obtained are more general 
than previous work because they can be used to describe strain-rate data for a range of 
temperatures. 
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The equations presented constitute a closed-form solution to the problem of iso- 
thermal secondary creep in thick-walled tubes loaded with either internal or external 
pressure. When the material properties, constants, and tube geometry are substituted 
into the equations, the creep strain rates and stresses are completely described 
throughout the wall of the tube. 

Fuel-pin cladding and heat-exchanger tubing are generally designed for thousands 
of hours of operation. The principal source of long-term deformation of these tubular 
structures will be secondary creep from pressures acting on the structures. 

Some observations can be made by comparing the equations for stress and strain 
rate for external and internal pressure loading. The equations for strain rate are 
similar but of opposite sign. The equations for stresses a , a Q , and cr are different, 
while the equations for the equivalent stresses are the same but of opposite sign. If the 
internal pressure and the external pressure acting on a tube are equal, the net equiva- 
lent stress will be zero everywhere in the tube wall. Consequently, no creep will occur. 


CONCLUDING REAAARKS 

Temperature-dependent strain-rate equations are presented for thick-walled tubes 
undergoing isothermal secondary creep from internal pressure. Temperature- 
dependent stress and strain-rate equations are provided for thick- walled tubes under- 
going isothermal secondary creep from external pressure. 

These equations are applicable to the analysis of stresses and long-term isothermal 
creep deformation of thick-walled tubular structures loaded by static internal or ex- 
ternal pressures. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 10, 1971, 

126-15. 
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APPENDIX - DERIVATION OF STRESS AND STRAIN-RATE EQUATIONS FOR 
THICK-WALLED TUBES AT CONSTANT TEMPERATURE 
UNDERGOING SECONDARY CREEP 
FROM INTERNAL PRESSURE 

The derivations in this appendix are from the analytical precedure presented by 
Bailey in reference 1. The assumption that = 0 is in agreement with experimental 
evidence. The assumption implies that planes perpendicular to the axis of a tube re- 
main plane. This is the condition of plane strain. 

Force equilibrium for cylindrical geometry in the absence of body force yields 


= a 0 - °r 


The strain displacement relations are 


r 


Differentiating with respect to time yields 


(where dots indicate d/dt). The compatibility equations for this case reduct to 


~ e r ~ e 9 
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Creep is a constant volume process, so that 


e r + e 0 + e z = 0 

Applying the assumption that e z = 0, we have 

e r = ~ e 6 

Substituting equation (A4) into (A3) yields 

d£ g _ 2 dr 

e 0 


(A4) 


Upon integration, we have 

k Q = Cr' 2 (A5) 

where the constant C is to be evaluated. 

The assumption that the principal shear strain rates are proportional to the princi- 
pal shear stresses can be expressed as follows: 



Substituting equation (A4) into (A6) yields 


o 


z 


a 9 + g r 
2 


(A6) 


(A7) 


An expression for the equivalent stress a is given by the distortion energy theory 
(also known as the von Mises yield criterion): 


a 





+ < a z 



(A8) 
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Substituting for a from equation (A7) gives 

Zi 


a " " a r ) 


The equivalent strain rate e corresponding to the equivalent stress is given by 

1 1/2 


6 E € r ' e 0 )2 + " e z )2 + (e z ' e r )2 ] 


Substituting for k z - 0 and e r = -e^ reduces equation (A10) to 

e = v !' 0 

Make the assumption that experimental strain-rate data can be represented by 

? = Bo«e- iH / RT 


Then 



= B 


& 

2 





n e -AH/RT 


and from equations (A5) and (Al) 



da 


/ 2 \(n+l)/n 



AH/RT 


Br 


1/n 


dr 

r 


Let 


D = 



Ce AH/RT 


B 


1/n 


(A9) 


(A10) 


(All) 


(A12) 


(A12) 


(A13) 
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da r = D 


dr 

r (n+2)/n 


Integration yields 



— + F 
2/ r 2/n 


The constants of integration C and F are evaluated using the boundary conditions . 
At r = b, a r =0 



_J_ + F = 0 
2/ b 2/n 


Dn 

2b 2 / n 


Dn 


2 1 

\b 2 / n 



At r = a, cr r = -p 


/Dn\ 

1 1 

\ 2 / 

b 2 / n a 2 / n _ 


With D replaced 




/ 2 \ n+1 Ce AH/RT ~ 

1/n (> i\ 

M B J 

\b 2/n a 2/ V 


-p 


Then 


il/n _ 


n 


P 


2 \ n+ - e 


1 o AH/RTl 1//n 



Substituting for C and letting x = b/r and p - b/a give 




2/n 


2/n _ 


(A15) 


Differentiating with respect to r in equation (A15) gives 


dg r MM x 2/n P 


dr 


W W p 2/n . j 


Substitution into the equilibrium equation (Al) yields the equation for : 

/— f>\ o /„ 

1 


cr„ = \ ■ P ^ p 


,2/n 


(A16) 


- 1 


From equation (A7) 


<r Z = 2 < 0 r + l 'e> 




1 - ELzJ. x 2/n 




n 


2/n 


- 1 


From equation (A9) 


x 2/n 

n p2/ n - 1 


p 


(A17) 


(A18) 
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From equations (A4), (All), and (A12) 


= -jt3 B5 n e -AH/BT 
r 2 

(A19) 

k„ =H B a"e- AH / RT 
Q 2 

(A20) 

, = 0 (assumption (5)) 

(A21) 


(A22) 
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